Digital three-state adiabatic passage 
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We explore protocols for three-state adiabatic passage where the tunnel matrix elements are varied 
digitally, rather than smoothly as is the case with conventional adiabatic passage. In particular, we 
focus on the STIRAP and related three-state schemes where the control is applied stepwise, with 
either equal spaced levels for the tunnel matrix elements or uniform pulse lengths. Our results show 
that the evolution typically shows the hallmarks of conventional adiabatic passage, although with 
additional resonances exhibiting no state transfer. 

PACS numbers: 05.60.Gg, 42.50.Ex, 03.67.-a, 03.65.Aa 
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I. INTRODUCTION 

Adiabatic passage is an interesting and important tool for 
manipulating quantum states. With adiabatic processes, 
one attempts to find a control strategy that maintains 
the system in an instantaneous eigenstate, whilst that 
Hamiltonian is smoothly varied from some initial state 
to the desired final state. Such control strategies have 
been used and proposed for a range of tasks including the 
manipulation of population within an atomic system [l|, 
0] , control of chemical reactions Q , control of vibrational 
population and adiabatic quantum computation Q. 

One of the tasks that adiabatic passage is particularly 
suited for is state transfer. State transfer is an important 
task in its own right, but also serves as a useful metric for 
evaluating new techniques. Here we focus explicitly on 
the task of state transfer in three-state systems using the 
STIRAP (STImulated Raman Adiabatic Passage) [J and 
related protocols such as Coherent Tunneling Adiabatic 
Passage (CTAP) [fHHl, Dark State Adiabatic Passage 
fDSAP 1 ) ll3l[T3 |. and adiabatic nonlinear frequency con- 
version [201, see Fig.[T] In these schemes we have an effec- 
tive three-state system where the states are denoted |1), 
1 2), 1 3) with nearest neighbour tunnel matrix elements 
Q,\ between states |1) and |2), and 0,2 between states |2) 
and 1 3). The task of these protocols is to effect popula- 
tion transfer from state |1) to |3) via the dark state, or its 
equivalent, \V ) = (0 2 |1) - O x |3>) /V^i + . Our aim 
here is to discover the limitations of employing piecewise 
constant, i.e. digital, tunnel matrix elements, rather than 
the smoothly varying control that is typically thought to 
be required. 

Piecewise adiabatic passage (PAP) was introduced by 
Shapiro et al. fl5| in the context of femtosecond pulse 
control of the internal states of atomic systems. In this 
case the control fields are applied for short periods of 
time with an overall envelope that ensures the evolution. 
Although formally this scheme cannot be considered as 
truly adiabatic due to the pulsed control, when viewed 
in the frequency domain, and observing that population 
cannot leak between eigenstates when the control fields 
are off, [ltjj . adiabatic-like evolution can be obtained. 
Controls comprising sequences of short pulses have been 
explored quite extensively recently in the context of PAP, 



see for example [l7|, [l8[ . The limits to such evolution in 
a general setting have also been explored by Boixo and 
Somma fl9l ]. 

Here we show the effects of employing a digital variation 
in the control fields, i.e. where the tunnel matrix ele- 
ments are forced to vary in discrete steps, rather than 
smoothly. This case is especially important for adiabatic 
frequency conversion using periodically poled media [2(| , 
and the inevitable digitisation that arises from the use 
of digital to analog converters for electronic control sig- 
nals. It is also likely that many computational models 
that use some form of finite clement method in solving 
adiabatic evolution show some de facto digitisation of 
the style we consider here. Our results show that digi- 
tal control can yield state transfer approximating more 
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FIG. 1: Three-state adiabatic passage has been explored in 
many systems with similar Hamiltonians, and four examples 
are shown here, (a) STIRAP effects state transfer between 
states of a three-level atom, typically ground states of a A 
system, where the tunnel matrix elements are the Rabi fre- 
quencies of control lasers, (b) A three-waveguide system can 
be used to effect CTAP variation where the tunnel matrix 
elements are varied due to the evanescent coupling between 
waveguides (c) CTAP in a triple-dot or triple-donor system 
is effected by gate-controlled variations in the wave function 
overlap of the particle between the sites, (d) For nonlinear fre- 
quency conversion, the required coupling variation is achieved 
by superimposing two chirped poling functions in a periodi- 
cally poled nonlinear material, assisted by two strong pump 
fields. 
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conventional three-state adiabatic passage, although we 
identify resonances where the transport does not exhibit 
adiabatic-like evolution. 

This paper is organised as follows. We first introduce the 
Hamiltonian and control scheme for the digital adiabatic 
passage (DAP). We then provide a simple analytic over- 
estimate of the fidelity of DAP. We consider two means 
for digitising the adiabatic passage, with uniform spac- 
ing of the pulses in time and uniform distribution of the 
tunnel matrix elements. Both of these cases are investi- 
gated, and a method for generating DAP schemes where 
there is uniform spacing in the control parameter is also 
discussed. Our results show adiabatic-like evolution with 
resonances of no evolution, and these resonances are used 
to identify optimal evolution strategies within the DAP 
framework. 



II. HAMILTONIAN 

We consider a three-state system described by a tight- 
binding model, where each state is coupled only to its 
neighbour. This model is useful for exploring STIRAP- 
like processes, although extensions to more realistic 
Hamiltonians have also been considered, especially in the 
context of spatial variations (i.e. CTAP), sec for exam- 
ple [ID, [2l|, [22[ . Although these adiabatic processes only 
require energy degeneracy between the end states, for 
simplicity we assume that all states have the same en- 
ergy. In this case the Hamiltonian becomes (with H = 1, 
and fii and real) 

W = n 1 (|l)<2| + |2)(l|)+n 2 (|2>(3| + |3)(2|). (1) 

Note that fli and are assumed to vary according to 
some time- varying external control, as in the case of STI- 
RAP and CTAP, or spatially varying design, as in the 
case of adiabatic passage in waveguides or nonlinear fre- 
quency conversion with periodically poled materials. For 
simplicity and without loss of generality, we will assume 
that the fl vary as a function of time. 
To understand the adiabatic evolution of the three-state 
system, we first solve for the eigenvectors of the Hamil- 
tonian, which are 



sin/cos 
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Wn) = 



± y/Qf + Ol|2) + Q 2 |3) 
flail) -ni|3) 



with corresponding eigenenergies 



E± 



fll = ±S, E = 0. 



(2a) 
(2b) 

(3) 



Here \T>o) is usually referred to as the dark state in the 
context of STIRAP, or the null state for CTAP. 
For adiabatic passage, the strategy to realise state trans- 
fer from |1) to 1 3) is relatively straightforward. One typi- 
cally chooses some smoothly-varying control function for 
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FIG. 2: Coupling schemes and eigenspectra for the case of 
sin/cos and linear coupling, (a) sin/cos coupling scheme 
showing both the continuous and digital variation with 15 
steps, (b) linear coupling scheme with showing both contin- 
uous and digital coupling for 15 steps, (c) Eigenspectrum for 
the sin/cos scheme. Note that the eigenvalues do not change 
throughout the protocol although the state compositions do 
change for each set of tunnel matrix elements, (d) Eigenspec- 
trum for the linear coupling scheme, showing that the point 
of minimum energy separation is at the midpoint of the pro- 
tocol. This is also the point where the adiabaticity parameter 
is maximised in the continuous version, (e) Adiabaticity for 
the sin/cos (green dashed line) and linear (red solid line) pro- 
tocols as a function of fractional time. Here we normalise the 
adiabaticity by multiplying by the total time. The sin/cos 
protocol is a straight line, whereas the linear protocol shows 
the more standard maximum in the adiabaticity at the mid- 
point of the evolution. 



fli and O2, with the constraints that H,2(t = 0) 3> f2i(i = 
0) and Qi(t = i max ) ^ = imax)i and the time t 

varies from to t max , where there is considerable overlap 
of the pulses for all intervening times. The functional 
form is relatively unimportant and many different pulse 
sequences have been e mployed . Here we employ a sinu- 
soidal control sequence [2j 



24| such that 



ill = VLm sin 



in 



2t„ 



f^2 = Am cos 



tlT 
— ' nia 



(4) 



as shown in Fig. [21a), where Am is the maximum tunnel 
matrix element. This control sequence has the advantage 
that none of the eigenvalues of the system vary during the 
adiabatic passage, i.e. E±(t) = ±Qm [Fig. Etc)]. 
The adiabaticity parameter is a measure for how well a 
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system stays in an instantaneous eigenstate as the Hamil- 
tonian is varied [25| . So for any two instantaneous eigen- 
states tpi and ip2 we define 



A = 



{MdtU\4> 2 ) 



KVi«i>-<V> 2 »2>p 



(5) 



and for adiabatic evolution, we require A <C f . For the 
pulse sequences typically employed for STIRAP/CTAP, 
the adiabaticity peaks either at or near the midpoint of 
the transfer [26[ (for a comparison of different schemes 
see [H|). For the case in cq. 01 however, the adiabaticity 
parameter is constant, namely 



As 



(6) 



We also consider a linear variation in the couplings 
[Fig. lib)], 

f^i = flM(t/t max ), fi 2 = — t/t max ). (7) 

This scheme has the advantage of simplicity and also cor- 
responds to systems where the control levels are equally 
spaced, but this scheme is less desirable for high fidelity 
transport, as will be shown below. The adiabaticity in 
the linear case is 



1 



t 



t 



1 



-3/2 



(8) 



In Fig. [5] we compare the control sequences, eigenener- 
gies and adiabaticity (in the smoothly varying case) for 
the two coupling schemes. Note that the linear evolu- 
tion shows the point of minimum energy separation at 
the midpoint of the protocol, in keeping with most of 
the other pulsing schemes that have been considered in 
the literature to date, whilst the eigenergies are constant 
for the sin/cos sequence. The fact that the adiabaticity 
for the linear sequence is larger for the same total time 
shows immediately that the linear scheme should be less 
efficient for state transfer than the sin/cos scheme. The 
constancy of the eigenenergies in the sin/cos scheme also 
has benefits for the DAP, as discussed below. 



III. EQUAL PULSE LENGTH DIGITISATION 

The aim of this work is to understand the consequences of 
digital control on three-state adiabatic passage protocols. 
It is important to realise that the concept of adiabaticity 
is formally inapplicable in this the assumption 

of ideal digital control implies step function variation in 
the control parameters. Hence adiabatic following is not 
strictly possible. Nevertheless, as we will show, the evo- 
lution obtained by a naive digitisation of the controlling 
tunnel matrix elements provides evolution that typically 



mimics the adiabatic evolution found for smoothly vary- 
ing control fields, as was the case with PAP fl5| . 
For digital control, the tunnel matrix elements are piece- 
wise constant with time. For N steps, we use 



ill = sin 



where 



2(N - 1) 



Qm cos 



2(N - 1) 



Nt 



(9) 



(10) 



and t G [0,i max ). 

If the evolution were adiabatic, then we could use adia- 
baticity arguments to determine the overlap of any given 
state with the eigenstate of the infinitcsimally incre- 
mented Hamiltonian. However, the digitisation intro- 
duces a series of discrete jumps. To determine an es- 
timate of the error, we therefore look at the overlap be- 
tween the null states before and after a change in the 
Hamiltonian, reasoning that in the limit of infinitely 
many steps that the continuously varying Hamiltonian 
should be approximated. Due to the simple form of the 
digital control, the difference between the £ th and (£+l) th 
null state is 



r] (Z) = l-\(V (0\V (Z + m 

IT 



2 (AT - 1) 



(11) 



In the limit that t|«l, the total error in the DAP is the 
sum over all the individual errors 



JV-l 

vt = vffi = N sil 



_2(N-1) 
which becomes in the limit of large N 



Vt = tt 2 /(4A). 



(12) 



(13) 



The above approximation gives a quick measure of the 
fidelity of the digital adiabatic process, confirming the 
intuition that increasing the number of steps should im- 
prove the overall fidelity. However, it misses some ex- 
tremely important physics, in particular the fact that 
population may be projected back into the null state 
from the other two eigenstates [27|, and also the fact 
that there will be interference in the accumulated phases 
of the residual population. Hence the error determined 
in eq. is in a real sense 'timeless'. The sin/cos evo- 
lution is particularly useful for understanding the evolu- 
tion, as all eigenstates have constant energy throughout 
the protocol, although of course they change their exact 
composition. 
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FIG. 3: Temporal evolution through DAP. Adiabatic-like evolution as a function of fractional time for (a) N — 5, t max = 5iv/Qm ', 
(b) N = 15, t max = 157r/f2jv/; and (c) N — 45, £ max = 457r/fiA/. The evolution is highly reminiscent of conventional three- 
state adiabatic passage, with the similarities becoming greater with increasing number of steps. When the time per step is 
r = 2nn/SlM then the evolution operator becomes the identity at the end of each step, hence the total evolution is the identity 
operator. This is illustrated for the cases that (d) N = 5, t m ax = Wtv/CIm', (e) N = 15, £ max = 307r/f2j\/; and (f) N = 45, 
f max = 907T /Q.M- Pseudocolor plots showing the evolution as a function of fractional time and total time for (g) N = 5, (h) 
N = 15 and (i) N = 45. Notice that the evolution is mostly adiabatic-like, with the addition of the resonances where the 
evolution is the identity. 



To perform accurate numerical calculations, we generate the unitary evolution operator corresponding to DAP, which 
is straightforward due to the piecewise constant nature of the digital Hamiltonian. The unitary evolution operator 
is constructed so that \ip(£t max /N)) = U(£)\ip[(£ — l)t max /N]). Note that for the sin/cos sequence the energy of all 
eigenstates is constant and so in this case £ = CIm- 



u(0 = i 



Ofcos(fr) i£Q,\ sin(£ t) O1O2 [cos(£r) — 1] 
i£Qi sin(£ r) £ 2 cos(£t) i£fl2 sm(£r) 
flifl 2 [cos(£ t) — 1] i£Q'2 sin(£ r) O 2 + 0| cos(£ r) 



(14) 



where we have introduced r, the length of the evolution under the current Hamiltonian. For equally spaced time steps 
t = t max /N. The total evolution over the DAP is 



U = U(N)U(N -!)••• [/(C) • • • U(2)U(1). 



(15) 



We first illustrate the process by showing several in- 
stances of the evolution for varying number of steps and 
total time. Fig. [3]shows the evolution of the populations, 
Pi = \(i\i)\ 2 , for TV = 5, 10, and 15. First considering the 
evolution through the protocol for varying £ max [Figs. [3] 
(g) (h) and (i)] we notice evolution that is strongly remi- 
niscent of adiabatic passage, with fairly smooth increase 
in population in |3). However there are also pronounced 



resonances where the final population in |3) is zero. The 
separation (in terms of i ma x) between these resonances 
increases as the number of digital steps is increased, and 
this result will be explained below. 

Figures [3] (a) (b) and (c) show examples of adiabatic- 
like evolution that is indicative of most of the parameter 
space. These results show that although the DAP does 
not show true adiabaticity, nevertheless the operational 
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FIG. 4: Final population in |3) after DAP using the sin/cos 
protocol as a function of t max for the first four non-trivial N 
for (a) N odd (N = 3, 5, 7, and 9) and (b) N even (N = 4, 6, 
8, and 10). When cos(Qmt) = 1, then the evolution becomes 
the identity, which is seen by the repeating nulls in the pop- 
ulation. Also of interest are the local maxima in population 
transfer, which correspond to the periodic Dirichlet kernel. 
As the number of digital steps increases, the number of local 
maxima increase and the pattern approaches a square func- 
tion. We have not shown the corresponding figures for the 
linear case as the results are less insightful. 

advantages of adiabatic passage (robustness against er- 
rors, smooth population variation and minimisation of 
population in the intervening state) are all preserved, 
with increasing number of steps improving the approxi- 
mation to true adiabatic processes. Notice that the max- 
imum population in \2) is given by rj(^) from Eg. 111! 
The resonances of low-fidelity transport are explored in 
Figs. |3] (d) (c) and (f), which show evolution through 
the centre of the resonances. Here we observe com- 
plex behaviour that does not transfer population to |3). 
We understand these resonances by noting that when 
tmax = InNir /Q,m , with n E N° in each case, then 
t = 2mr/Q,M and E/(£) is the identity matrix. Hence 
the overall evolution must also be the identity. 
In Fig. 2] we plot the population in state |3) at the end of 
the DAP protocol as a function of t max for different val- 
ues of N . The resonances corresponding to when hi = I 
are seen by the dips as a function of t max . Also of interest 
are the local maxima and minima in population transfer 
between these resonances. We find that the pattern of 
these extrema phenomenologically corresponds to a sub- 
set of the periodic Dirichlet kernel, however the heights 
of the nodes did not. We also note that protocols with 
odd TV and total time of Ntt, exhibit comparable evolu- 
tion to those reported in Ref. [llj]. The elegance of the 
results for the sin / cos scheme are not replicated with 
other protocols for the timing scheme utilised. 
We further explore the fidelity by showing in Fig. [5] the 
error in DAP transfer as a function of N and t max for 
both the sin/cos and linear schemes. The traces in Fig.0] 
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FIG. 5: Transport error [1 — |(^(imax)|3)| ] as a function of 
N and i ma x, for (a) the sin/cos scheme, and (b) the linear 
scheme. There is relatively low error away from the identity 
evolution resonances and with increasing N the spacing of 
these resonances increases linearly as does the number of op- 
timal nodes. In each case, the horizontal axis begins at N=3 
as this is the first non-trivial result. Note that the sin/cos 
scheme has sharper resonances as the eigenergies are con- 
stant throughout the protocol, whereas the changing energy 
observed in the linear scheme leads to a band over which some 
of the evolutions are the identity. 



correspond to vertical slices through Fig. [S] (a). Again, 
noticeable are the resonances showing identity evolution, 
and also apparent are the local maxima identified above. 
For the linear scheme [Fig[5](c)], the resonances are more 
complicated and do not result in a perfect identity. This 
is because 17 f + Vt\ is n0 longer constant, and hence per- 
fect identity evolution across every digital step is not pos- 
sible for a given r. 

We also compare the final state population for even and 
odd N with the analytical estimate from Eq. [T2] for in- 
creasing N in Fig. [5] The analytic estimate does not 
capture the resonant phenomena, but does show the in- 
creasing fidelity of DAP with increasing N. In both the 
linear and sin/cos schemes we see a general trend in the 
reduction of error with increasing N and a difference be- 
tween evolution where the total time is an even or odd 
multiple of tt/Qm- However, as was evident in Fig.[5l the 
poor fidelity resonances are less clear in the linear scheme 
than for the sin/cos scheme. This result is a consequence 
of the fact that the eigenspectrum has different values as 
a function of fractional time in the linear scheme, and 
hence the resonances are broader and less pronounced. 
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FIG. 6: Transport fidelity at t max = 100tv/Q,m (green) and 
fmax = Win/ SI m (blue) for varying number of steps, N, com- 
pared with the analytical results for n (red), (a) Shows results 
for the sin/cos scheme, and (b) shows the case for the linear 
scheme. As expected, as N increases, the fidelity improves, 
except for resonance points as discussed above. Notice the 
difference between the even and odd results, which is more 
pronounced for the sin/cos scheme than the linear scheme. 



IV. UNIFORMLY VARYING TUNNEL MATRIX 
ELEMENTS, NON-UNIFORM PULSE LENGTH 

In the previous section we considered the simplest ap- 
proach with respect to timing; the pulses defining the 
protocol were designed to have equal temporal extent. 
Although conceptually straightforward, in some situa- 
tions the control over the timing of pulses may be greater 
than the control over the level of the pulses, and hence 
for such situations it is more natural to consider dig- 
ital schemes where the pulse length is optimised for a 
given digitisation level. As the ultimate limit of such 
control imperatives, we may consider bang-bang control 
as a one-bit control scheme with fine timing control [28| ■ 
For simplicity we consider only the case where the tunnel 
matrix element is varied between uniformly spaced lev- 
els, however our method may be easily generalised to the 
case where the parameter controlling the tunnel matrix 
element (e.g. a voltage level for gate-driven processes or 
waveguide separation for waveguide adiabatic passage) is 
the uniformly spaced digital quantity. 
To design the appropriate pulse length, recall that for 
a given digital step, identity evolution occurs when the 
pulse length is r = 2n/£. Optimal complete transfer 
occurs for step £ when the length corresponds to 

rz = j, (16) 

where £ = £ (£). More generally, one can determine an 
equivalent to Fig. 0] for the compensated linear scheme, 
which is stretched by a factor dependent on the number 
of steps and the precise details of the scheme used. 




FIG. 7: (a) Linear scheme with energy-dependent pulse 
length. The tunnel matrix elements are taken from a set of 
equally spaced levels and the pulse time is inversely propor- 
tional to the energy, (b) Variation in pulse length relative 
to the uniformly spaced scheme. More time is spent in lower 
energy areas than high meaning that the pulse length at the 
midpoint of the protocol is longer than at any other time in 
the protocol. 



The compensation scheme described above provides a 
method to determine the pulse length for any control pro- 
tocol where the £ are known. To illustrate this method 
for setting the pulse length, we return to the linear con- 
trol scheme where wc require the tunnel matrix elements 
to be varied between equally spaced levels. However we 
now optimise the time according to eq. UHl This modified 
control scheme is illustrated in Fig. [7j which shows the 
control protocol as well as the shift in the pulse length rel- 
ative to the uniform pulse time protocol discussed above. 
The optimal pulse time approach means that the pulse 
length at the midpoint of the protocol is longer than at 
the beginning or end of the protocol. Interestingly, this 
is a similar result to that obtained from standard adia- 
baticity analysis, although the reason for longer pulses is 
different in this case. 

In Fig. [5] we compare the evolution through DAP for 
the linear scheme without correction to the pulse time 
(i.e. uniform pulse length) [for N = 7 in Fig. M, a ) and 
N = 45 in Fig. M, c )] and with the correction applied [for 
N = 7 in Fig. Mb) and N = 45 in Fig. MA)]- Popula- 
tions are shown on a logarithmic scale to highlight the 
transient population in |2). As expected, the evolution 
is far more regular for the temporally compensated cases 
than the uniform pulse time results. The regularisation 
of the evolution is also seen in Fig. |H1 which shows trans- 
port fidelity as a function of N and t max , and should be 
compared with the plots in Fig. [5j Note that the reso- 
nances of identity evolution have sharpened so as to be 
comparable to those of the sin/cos protocol. 



V. CONCLUSIONS 

Digital adiabatic passage, when applied to three state 
systems, is shown to approximate conventional adiabatic 
passage in the limit of large N. Although this was ex- 
pected, it is perhaps surprising that the adiabatic limit 
should be well approximated even for relatively few digi- 
tal steps, albeit with certain restrictions on the total time 
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FIG. 8: Population through the linear protocol showing the 
effects of temporal compensation, (a) Uncompensated evolu- 
tion with N = 7, (b) compensated evolution with N = 7, 
(c) uncompensated evolution for N = 45, and (d) com- 
pensated evolution for N = 45. In each case we have set 
imax = T (£)> i- e - we have the set the total time to that of 
the compensated protocol. Both (a,c) have residual popula- 
tions in |1) and 2) at t max , and also erratic transient popu- 
lation in |2) whereas (b,d) have complete population transfer 
with more regular transient |2) population. 




Number of steps, N 

FIG. 9: Transport fidelity for the linear scheme with com- 
pensated pulse times. This figure shows that the compen- 
sated pulses smoothens the function and gives rise to superior 
transport than the equally spaced protocol [c.f. Fig. [5jb)]. 



of the protocol. In addition to adiabatic-likc behaviour, 
we also find 'resonances' exhibiting poor transport, su- 
perimposed on the adiabatic-likc evolution. These reso- 
nances of poor transport arise due to evolution over some 
or all of the digital steps resulting in the identity. 
We compared two schemes for effecting digital adia- 
batic passage: sin/cos and linear variation; for uniform 
pulse lengths and compensated pulse lengths. Without 
pulse length compensation, the sin/cos scheme gave su- 
perior transport results. The overall transport was gen- 
erally higher fidelity as the adiabaticity parameter of the 
sin/cos scheme was less than that for the linear scheme, 
however also the resonances of poor transport were more 
clearly defined due to the fact that the eigenergies were 
constant throughout the protocol. Our results show that 
the constant eigenenergies of the sin/cos scheme provide 
improvements to the more commonly investigated adia- 
batic passage schemes. 

Optimisation of the pulse length for the linear scheme im- 
proved the fidelity of this approach to that of the sin/cos 
approach. Whilst non-uniform pulse lengths are more 
complicated, in some systems, it may be possible to em- 
ploy finer control over timing than the tunnel matrix ele- 
ment, and hence improve the fidelity of digital adiabatic 
passage. We assumed that the tunnel matrix elements 
were selected from a set of uniformly spaced levels. How- 
ever, our method may be generalised to cases where the 
tunnel matrix elements vary within any set of discrete 
levels. 

Overall, our results show that adiabatic passage is ro- 
bust against digitisation of the control parameters, a re- 
sult that is not obvious from simple adiabatic analysis. 
The implications of these results is that the class of sys- 
tems that are amenable to adiabatic passage techniques 
is larger than previously thought. 
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